Spontaneous fourfold-symmetry breaking driven by electron-lattice coupling and 

strong correlations in high-T c cuprates 
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Using dynamical-mean-field theory for clusters, we study the two-dimensional Hubbard model in 
which electrons are coupled with the orthorhombic lattice distortions through the modulation in the 
hopping matrix. Instability towards spontaneous symmetry breaking from a tetragonal symmet- 
ric phase to an orthorhombic distorted phase is examined as a function of doping and interaction 
strength. A very strong instability is found in the underdoped pseudogap regime when the inter- 
action strength is large enough to yield the Mott insulating phase at half filling. The symmetry 
breaking accompanies the recovery of quasiparticle weights along one of the two antinodal directions, 
leading to the characteristic Fermi arc reconnection. We discuss the implications of our results to 
the fourfold symmetry breaking reported in systems where the underlying crystal does not have any 
structural anisotropy. 

PACS numbers: 71.10.Fd,74.25.Jb,74.72.Kf 



Electronic ncmaticity has become one of the central 
subjects of correlated-clcctron systems^ - — For high-T c 
cuprates, very large anisotropics in low energy excitations 
have been experimentally reported and their connec- 
tion with the "pseudogap phase" has been discussed. 

In a system such as Yl^CusOg+a; (YBCO), there 
exists an intrinsic structural anisotropy resulting in a 
tiny but finite band-structure anisotropy. This band 
anisotropy has been shown to induce huge effects in 
low-energy excitations^ when the system is close to a 
correlation-induced Pomcranchuk instability— Although 
cluster dynamical mean-field studies do not find a spon- 
taneous symmetry breaking in a two-dimensional single- 
band Hubbard model ; 15 ' 16 a tiny band anisotropy was 
shown to dramatically amplify the anisotropy in the 
dc transport and electronic excitation spectrum in the 
underdoped pseudogap regimei 16 ' 17 On the other hand, 
in a system such as Bi 2 Sr2Ca„_iCu„08+; / (BSCCO) 
and C^-zNa^CuC^C^ (CNCOC), there is no intrinsic 
structural anisotropy but the symmetry breaking from 
tetragonal (C4) to orthorhombic (C2), called intra-unit- 
cell (IUC) nematicity, has been observed] 11 ' 13 Within a 
mean-field treatment, IUC order in the Emery model for 
the Cu02 plane has been analyzed^ However, very large 
interactions arc required to realize the IUC symmetry 
breaking. This may indicate the importance of additional 
degrees of freedom. 

Here, we consider a correlated model for cuprates in- 
cluding the coupling between electrons and lattice distor- 
tions (EL) as a possible ingredient for experimentally re- 
ported spontaneous C4 symmetry breaking. We observed 
a moderate tendency towards symmetry breaking when 
the chemical potential is located near the van Hove singu- 
larity. In addition, we found a very strong instability in 
the underdoped pseudogap regime when the interaction 
strength is large enough to yield a Mott insulating state 
at half filling. The stabilization of the distorted phase 
comes from the gain in the kinetic energy In the over- 
doped regime, the Fermi surface is deformed to split the 
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FIG. 1: Model electron-lattice coupling. Orthorhombic dis- 
tortion u induces the hopping anisotropy with the coupling 
constant a. 



van Hove singularity at (it, 0) and (0, it) and shift it from 
the Fermi level, while in the underdoped regime in the 
presence of strong Coulomb interaction the pseudogap 
becomes anisotropic. Our results may provide a coherent 
picture connecting the "electron nematicity" and pseu- 
dogap behavior in high-T c cuprates highlighting the dif- 
ference between YBCO, whose structure is intrinsically 
anisotropic, and BSCCO, which is structurally isotropic 
but the C4 symmetry is found to be locally broken. 
Wc consider the following electron-lattice coupled 



model: H = H e j e + Hi at t- H e i e is the two-dimensional 
(2D) Hubbard model as a generic model for high-T c 
cuprates: 
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Here, di a is the annihilation operator for an electron 
with spin a at site i and U is the local Coulomb in- 
teraction. The band structure part is described by t^; 
^^■corresponds to the transfer integral and and U—j = fx 
is the chemical potential. The lattice part is given by 
Hiatt = -^-u 2 , where u is the orthorhombic distortion, 
K is the elastic constant, and N is the total number of 
sites. We consider the coupling between the orthorhom- 
bic distortion and electrons as the modulation in the 
nearest-neighbor (NN) transfer integral t along the x and 
y directions as t x ^ y = t ± au as illustrated in Fig. [TJ with 
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the next-nearest-neighbor transfer integral t' unaffected 
by the distortion. This electron-lattice coupling may also 
be regarded as a simplified one realized in, for example, 
a low-temperature tetragonal (LTT) phase. Here, the 
crystal structure is "tetragonal" with the equal lattice 
constants along the a and b axes but the electronic band 
structure is "orthorhombic" due to the coherent rotation 
of Cu0 6 octahedrai^^ 

We analyze our model using the cellular dynamical- 
mean-field theory (CDMFT) (Refs. HH and [H) at zero 
temperature. This method captures the full dynam- 
ics [i.e., the frequency dependence of the spectral func- 
tion (SF)] and the short-ranged spatial correlations be- 
yond the single-site dynamical-mcan-field theory and has 
been applied for a variety of problems in low-dimensional 
systems ,22r-2£ The CDMFT maps the bulk lattice prob- 
lem onto an effective Anderson model describing a clus- 
ter embedded in a bath of noninteracting electrons. The 
short-ranged dynamical correlations are treated exactly 
within the cluster. In this study, we employ a 2 x 2 pla- 
quette (N c = 4) coupled to eight bath orbitals and solve 
it using the Lanczos exact diagonalization technique^ 
which requires a low-energy cutoff w c corresponding to 
the discrete imaginary frequency as u) n = (2n + l)u> c . In 
this work, we take w c = 2x 10~ 2 7ri. ui c should not be in- 
terpreted as real temperature (times tt) because only the 
ground state of the impurity model is taken to compute 
the Green's function. The numerical details are described 
in Ref . [H 

In the following discussion, we use as a band parameter 
t' = — Q.3t that is appropriate for cuprates. For interac- 
tion strength, we consider U = lOi, At and 0. The largest 
U is supposed to be relevant for cuprates. 

We start from the discussion on the first instability 
caused by the EL coupling. This could in principle be 
done by computing the C2 susceptibility and finding a 
parameter range where the susceptibility diverges. Such 
an analysis normally requires the inclusion of vertex cor- 
rections, but the precise form is unknown for the current 
CDMFT technique. Instead, we utilize the Ginzburg- 
Landau theory with the self-energy functional approach 2 -^ 
by which the electronic contribution F e i e to the free en- 
ergy F t ot = F ete + Fi a tt is written as 

F ele = -T]T / lndetjl - [t(k) - t c - t(iLo n )]G c {iu n )} 
+F C . (2) 



„jj is the hopping matrix on the 2x2 plaquette, 



Here /'.. = / 

and /,,ik: = N^^e 



e K+k describes the 



hopping between the clusters covering the original lattice 
k are wave vectors in the reduced Brillouin zone, and 
K = (0, 0), (tt, 0), (0, 7r), and (n, tt). The bare dispersion 
is given by £k = — 2{t x cos k x +t y cos k y +2t' cos k x cos k y ). 
T{iui n ) is the hybridization function with which the clus- 
ter Green's function at the ground state is written as 
G c (iu> n ) = [iuj n ~t c — T(ioj n ) — E(iu; n )] , where T,(iu> n ) 
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FIG. 2: (Color online) 5 t dependence of electronic free energy 
AF e i e (a) and (b) and the double occupancy D (c) and (d). 
V = At for (a) and (c), and U = KM for (b) and (d). Solid 
lines in (a) and (b) are quadratic fits to the numerical data 
of AF ele . 



is the cluster self-energy. Because of the low-energy cut- 
off uj c , the free energy is approximately calculated using 
the finite-temperature form Eq. ^ at T — uj c /^- As a 
result, the results are somewhat smeared out, underesti- 
mating the instability. Finally, F c is the free energy of 
the cluster model. The lattice contribution is given by 



F 



latt 



In practice, we compute F e i e as a function of the 
hopping anisotropy S t = au as AF e i e (S t ) = F e i e (8 t ) 



F e i e (0) 



/35? + 0{5f). Normally, the 0(S?) contri- 



bution is positive. As the lattice contribution F[ att — > 
7^2 KS 2 is always positive, quadratic fitting to AF e i e (S t ) 
gives the critical EL coupling towards the spontaneous C4 
symmetry breaking as a 2 / K = 1//3. The critical point 
deduced in this way signals the second-order transition. 
When the 0(Sf ) contribution is negative, the transition 
becomes first order, thus, a 2 /K = 1//3 should be re- 
garded as the upper limit of the coupling above which 
a distortionless state is no longer an energy minimum. 
Instead, the critical EL coupling for the first-order tran- 
sition is located at the smallest a 2 /K which satisfies 
AF e ie{$t) + ^kK5? = at S t 0. Because of the higher- 
order terms 0(Sf) and 0(8f) in F e i e and the small change 
in the carrier density with S t , these fitting procedures are 
notoriously difficult. Nevertheless, an overall trend can 
be deduced. 

As examples, we plot AF e i e (St) in Figs.[2ja) and 2(b). 
AF e i e decreases with increasing St- According to the lin- 
ear combination of atomic orbitals method^^ the EL 
coupling constant a is estimated to be 0.9 eV/A, and the 
distortion considered here is rather small; St = 0.03 cor- 
responds to u ~ 0.033 A. In most cases, AF e i e deviates 
from the quadratic curve upwards. While for U = At 
with fj, = 0.8t and U = lOt with \x = 1.8t, AF e i e shows 
a small downward deviation, indicating first order tran- 
sitions. For U = lOt with ji = 2.0t, corresponding to the 
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FIG. 3: (Color online) Phase diagram for the 2D Hubbard 
model with the EL coupling as a function of electron density 
N and the coupling constant a 2 /K. Parameters are U = Wt 
(squares) and U = 4t (circles) with t' = — 0.3t. The critical 
points for the second- (first-) order transition are indicated by 
filled (open) symbols. The critical points where distortionless 
states lose metastability are indicated by crossed symbols. 
Thick lines are guides to the eye. For comparison, the phase 
boundary for U = is shown as a light solid line. 



undcrdopcd regime. A.F e i e shows a strong dependence 
on St, indicating strong instability. 

As shown in Fig. [3J the AF e i e -8t curve is rather sen- 
sitive to the doping concentration. To sec its origin, 
we plot, in Figs. [U[c) and 2(d), the double occupancy 
D = (d^difdl^dii) as a function of 6 t . One notices that 
the double occupancy remains unchanged, i.e., the po- 
tential energy UD remains unchanged. Thus, the gain in 
the "kinetic energy" dominates the AF e i e -6t irrespective 
of doping dependence. 

Figure [3] shows the resulting phase diagram^ - As de- 
scribed below, instabilities appear at two doping regimes 
due to different mechanisms, say type A and type B. 
Type A is a weak-coupling mechanism and appears for 
both U = At (boundary is indicated by circles) and 
U = lOt (squares) at N ~ 0.8, near the van Hove fill- 
ing. Type B, on the other hand, is a strong-coupling 
mechanism and only appears when U is large in the un- 
derdoped regime (N > 0.9 for U = lOt). 

Type-A weak coupling instability also appears for U = 
0. In this case, the critical coupling is given by expand- 
ing the free energy up to 0(6%) and equating its coef- 
ficient to zero, and the resulting expression is a 2 / K = 

— 7r 2 /2 fdk 2 (cosk x — cos/c y ) 2 /'(ek — n), where /' is the 
derivative of the Fermi-Dirac distribution function. Due 
to the logarithmic divergence in the DOS, the critical 
coupling is minimized at the van Hove filling N ~ 0.726 
as shown as a light solid line. By finite U, the instability 
is shifted to larger N ~ 0.8. From the analysis of SFs, 
the shift in the critical N is caused by the enhanced band 
anisotropy due to correlations. The corresponding SFs 
[contour plot of the SF, A(k,u = 0) = -±ImG(k, oj = 0)] 
are presented in Fig.|3](top left, isotropic band) and (mid- 
dle left, anisotropic band), where the Green's function 
pcriodization scheme is adopted^ The enhanced band 



anisotropy is evident from the comparison with the FS 
for U = (a white line). By finite U, the FS opens up 
to become quasi-one-dimcnsional [sec Fig. [4] (middle left 
panel)]. This is favorable for gaining the kinetic energy 
by splitting the van Hove singularity and shifting it from 
the Fermi level. This could also explain why the first- 
order transition appears at N > 0.8, where correlation 
effects are stronger <2i 

For the type-B instability, the band anisotropy re- 
mains almost unchanged but the anisotropy in the scat- 
tering rate is enhanced significantly. These points can be 
clearly seen in the corresponding SFs presented in Fig. @] 
(top right, isotropic band, and middle right, anisotropic 
band). Here more importantly, the coherence is recov- 
ered near (0, tt) by the band anisotropy because the FS 
goes away from the so-called "hot spot," while the FS 
near (tt, 0) approaches the hot spot. This results in the 
reconncction of the "Fermi arc" between the first and 
the second Brillouin zones neighboring along the y direc- 
tion; see Fig. [4] (middle right) and Fig. 1 (d) in Ref. [HI 
The relation between the quasiparticle coherency and the 
kinetic energy can be directly seen from the expression 
for the kinetic energy E kme = TJ2u. n £ kG(k, iu) n ). As 
the quasiparticle coherency is lost in the "symmetric" un- 
dcrdopcd pseudogap regime, the gain in the quasiparticle 
weight leads the dramatic gain in the kinetic energy as 
seen in Fig. H](b), resulting in the strong instability in 
this regime. 

Given the above discussion, the typc-B instability is 
expected to be more relevant for undcrdopcd high-T c 
cupratcs compared with the type-A instability and other 
weak-coupling instabilities. The type A instability is sup- 
pressed by correlations because the van Hove singularity 
is smeared out by the imaginary part of the self-energy. 
In fact, for U = lOt, this instability is almost diminished. 
Not only by the imaginary part of the self-energy, the 
type-A instability is also suppressed by finite tempera- 
ture. Because of the linear T dependence of the electronic 
free energy, the instability is expected to go away rather 
quickly at elevated temperatures. On the other hand 
the type B instability requires large U, resulting in the 
pseudogap or the Fermi arc which sets its energy scale. 
Therefore, the type B instability is expected to survive at 
relatively high temperatures as long as the pseudogap or 
Fermi arc remains. For U = lOt, there appears a dip in 
the instability at N ~ 0.88 because both instabilities are 
weak (see Fig. [3} at this doping. This does not contra- 
dict the seemingly stronger instability near 1/8 doping 
reported experimentally^ as it has different origins. Ad- 
ditional weak coupling instability incompatible with our 
model could take place at carrier densities smaller than 
N ~ 0.9 as discussed, for example, in Refs. [H andl37l 

When the system is in the vicinity of the structural 
transition, dynamical or statistical fluctuation effects 
should play important roles. If the transition is of the 
first order, such a critical regime would be characterized 
by a superposition between two distortion modes which 
minimize the free energy. Furthermore, even in a dis- 
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FIG. 4: (Color online) Evolution of the spectral function at 
the Fermi level in the first quadrant of the Brillouin zone 
for U — Wt. Left panels: results for N = 0.83 (overdoped 
regime), and right panels: results for N = 0.94 (underdoped 
regime). Top: Symmetric SF; middle: asymmetric SF with 
St/t = 0.03; bottom: averaged SF with S t /t = ±0.03. The 
Green's-function periodization scheme is used with the small 
imaginary part ir\ (r\ = O.li) added to the real frequency. The 
maximum of the spectral weight is A max t = 0.6 (left) and 0.35 
(right). A white line in the middle-left panel shows the FS 
for 5 t /t = 0.03, N = 0.83 and [7 = 0. 



torted phase, a sample could form domains. As a result, 
low-spatial-resolution angle-resolved photoemission spec- 
troscopy (ARPES) measurements would detect the SF 
that is an average of SFs over a finite lattice spacing. We 
simulate the latter two cases by taking the average of the 
SFs with different anisotropy parameters, St/t = ±0.03. 
Figure 2] (bottom left and bottom right) show the av- 
eraged SF for N = 0.83 (the overdoped regime) and 
N = 0.94 (the underdoped regime), respectively. For 
the overdoped regime, the spectral function is broadened 
at (7T, 0) and (0,7r) relative to the results without dis- 
tortion. Thus, the "lattice fluctuation" and structural 
domain formation act as if enhancing the pseudogap be- 
havior. A similar effect has been reported for the thermal 
nematic fluctuation^ On the other hand, for the under- 
doped regime, the SF in the undistorted phase and the 
averaged one in the distorted phase arc nearly identical. 
This is because the shape of the FS is insensitive to the 
band anisotropy in this doping regime. 



In contrast to YBCO, BSCCO, La 2 - a: Sr :c Cu04 and 
CNCOC do not have a source for band anisotropy, yet 
(local) C4 symmetry breaking has been reported. For 
BSCCO, recent scanning micro-x-ray-diffraction^ and 
scanning tunneling microscopy 4 ^ measurements revealed 
that the system is structurally inhomogencous involv- 
ing LTT-like distorted regions. These results are con- 
sistent with our picture if the EL coupling is in the 
range of the spontaneous distortion. In fact, from the 
measured elastic constants c ~ 1.7 x 10 12 dyn/cm 2 for 
La 2 Cu0 4 (Refjp) and c - 1.3 x 10 12 dyn/cm 2 for 
BSCCO (Ref. |42I ), our elastic constant is estimated as 
K ~ 6 eV/A 2 for both systems. The resulting EL cou- 
pling constant^ 9 . a 2 /Kt <~ 0.2 locates these systems in- 
side the spontaneous distortion regime. Thus, it is de- 
sirable to experimentally clarify the relation between the 
electroni o 11 ' 13 and the structural C4 symmetry breaking. 

Our prediction can be tested by high-spatial-resolution 
ARPES measurements as in the x-ray-diffraction mea- 
surements in Ref. [391 For underdoped cuprates below 
the pseudogap temperature, we expect that anisotropic 
FSs as in the lightly doped YBCO (Ref. |H) or recon- 
nected Fermi arcs are spatially distributed. The pseu- 
dogap behavior should preempt or accompany the local 
lattice distortion with small effects on the "bulk" SF at 
the structural transition. In contrast, the opening of a 
pseudogap and the local lattice distortion are expected to 
take place simultaneously in the overdoped regime. From 
these experimental tests, a variety of anomalies in rela- 
tion to the electronic nematicity and pseudogap behav- 
ior in high-T c cuprates could be coherently understood 
in terms of the EL coupling and the absence/presence 
of the intrinsic band anisotropy. Further, Raman scat- 
tering would be a useful tool to distinguish different 
roles played by electronic systems and lattice (phononic) 
system o 44 ' 45 because the electronic contributions to the 
Big Raman scattering intensity are suppressed in the 
pseudogap regime^ 

In addition to high-T c cuprates, electronic nematic- 
ity was suggested for the double-layer ruthenates un- 
der an applied magnetic field 5 -^ and Fe-based high-T c 
superconductors above magnetic phase transitions £^ 
A number of theoretical scenarios have been proposed, 
including microscopic phase separatio n 48 ' 49 and quasi- 
one-dimensional orbital ordering 5 ^ for ruthenates, and 
fluctuating magnetic stripe order ; 51 ' 52 a ferro-orbital 
ordering j 5 - 3 - and magnetoelastic couplin g 54 ' 55 for Fe su- 
perconductors. It is worth investigating these scenarios 
including EL coupling with electron correlations treated 
beyond static mean-field approximations. 

Summarizing, we investigated the C4 symmetry break- 
ing by the coupling between electrons and the lattice dis- 
tortion using an interacting model for cuprates within 
CDMFT. We found the strong instability towards C4 
symmetry breaking in the underdoped pseudogap regime 
in the presence of the strong interaction yielding the Mott 
transition. Thus, this instability is a strong-coupling ef- 
fect characteristic of a doped Mott insulator. Addition- 
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ally, a weak-coupling instability exists near the van Hove 
filling, but this instability is suppressed by strong corre- 
lations because the imaginary part of the electron self- 
energy smears the van Hove singularity. On the other 
hand, the imaginary part of the self-energy plays an es- 
sential role for the strong-coupling symmetry breaking 
by increasing the kinetic energy by recovering the quasi- 
particle coherence. This leads to the characteristic Fermi 



arc reconnection below the pseudogap temperature. Our 
finding can be tested by high-spatial-resolution ARPES. 

The authors acknowledge discussions with V. R. 
Cooper, H.-Y. Kee, and H. Yamase, and the support from 
the U.S. Department of Energy, Basic Energy Sciences, 
Materials Sciences and Engineering Division (S.O.), and 
the Grants-in-Aid for Scientific research, MEXT, Japan 
(N. F.). 



* okapon@ornl.gov 

1 H. Yamase and H. Kohno, J. Phys. Soc. Jpn. 69, 332 
(2000); 69, 2151 (2000). 

2 C. J. Halboth and W. Metzner, Phys. Rev. Lett. 85, 5162 
(2000). 

3 S. A. Kivelson, I. P. Bindloss, E. Fradkin, V. Oganesyan, J. 
M. Tranquada, A. Kapitulnik, and C. Howald, Rev. Mod. 
Phys. 75, 1201 (2003). 

4 E. Fradkin, S. Kivelson, M. Lawler, J. Eisenstein, and A. 
Mackenzie, Annual Review of Condensed Matter Physics 
1, 153 (2010). 

5 R. A. Borzi, S. A. Grigera, J. Farrell, R. S. Perry, S. J. 
S. Lister, S. L. Lee, D. A. Tennant, Y. Maeno, and A. P. 
Mackenziel, Science 315, 214 (2007). 

6 J.-H. Chu, J. G. Analytis, K. De Greve, P. L. McMahon, 
Z. Islam, Y. Yamamoto, I. R. Fisher, Science 329, 824 

(2010) . 

7 Y. Ando, K. Segawa, S. Komiya, A. N. Lavrov, Phys. Rev. 
Lett. 88, 137005 (2002). 

8 Y.-S. Lee, K. Segawa, Y. Ando, D. N. Basov, Phys. Rev. 
B 70, 014518 (2004). 

9 V. Hinkov, D. Haug, B. Fauque, P. Bourges, Y. Sidis, A. 
Ivanov, C. Bernhard, C. T. Lin, B. Keimer, Science 319, 
597 (2008). 

10 R. Daou, J. Chang, D. LeBoeuf, O. Cyr-Choiniere, F. Lal- 
iberte, N. Doiron-Leyraud, B. J. Ramshaw, R. Liang, D. 
A. Bonn, W. N. Hardy, L. Taillefer, Nature (London) 463, 
519 (2010). 

11 M. J. Lawler, K. Fujita, Jhinhwan Lee, A. R. Schmidt, Y. 
Kohsaka, C. K. Kim, H. Eisaki, S. Uchida, J. C. Davis, J. 
P. Sethna, E.-A. Kim, Nature (London) 466, 347, (2010). 

12 D. Haug, V. Hinkov, Y. Sidis, P. Bourges, N. B. Chris- 
tensen, A. Ivanov, T. Keller, C. T. Lin and B. Keimer, 
New J. Phys. 12, 105006 (2010). 

13 Y. Kohsaka, T. Hanaguri, M. Azuma, M. Takano, J. C. 
Davis, and H. Takagi, Nature Phys. 8, 534 (2012). 

14 H. Yamase and W. Metzner, Phys. Rev. B 73, 214517 
(2006). 

15 E. Gull, O. Parcollet, P. Werner, and A. J. Millis, Phys. 
Rev. B 80, 245102 (2009). 

16 S. Okamoto, D. Senechal, M. Civelli, A.-M. S. Tremblay, 
Phys. Rev. B 82, 180511(R) (2010). 

17 S.-Q. Su and T. A. Maier, Phys. Rev. B 84, 220506 (2011). 

18 M. H. Fischer and E.-A. Kim, Phys. Rev. B 84, 144502 

(2011) . 

19 J. D. Axe, A. H. Moudden, D. Hohlwein, D. E. Cox, K. 
M. Mohanty, A. R. Moodenbaugh, and Y. Xu, Phys. Rev. 
Lett. 62, 2751 (1989). 

20 B. Nachumi, Y. Fudamoto, A. Keren, K. M. Kojima, M. 
Larkin, G. M. Luke, J. Merrin, O. Tchernyshyov, Y. J. 
Uemura, N. Ichikawa, M. Goto, H. Takagi, S. Uchida, M. 



K. Crawford, E. M. McCarron, D. E. MacLaughlin, and 
R. H. Heffner, Phys. Rev. B 58, 8760 (1998). 

21 G. Kotliar, S. Y. Savrasov, G. Palsson, and G. Biroli, Phys. 
Rev. Lett. 87, 186401 (2001). 

22 G. Kotliar, S. Y. Savrasov, K. Haule, V. S. Oudovenko, O. 
Parcollet, and C. A. Marianetti, Rev. Mod. Phys. 78, 865 
(2006). 

23 M. Capone, M. Civelli, S. S. Kancharla, C. Castellani, and 
G. Kotliar, Phys. Rev. B 69, 195105 (2004). 

24 B. Kyung, S. S. Kancharla, D. Senechal, A.-M. S. Trem- 
blay, M. Civelli, and G. Kotliar, Phys. Rev. B 73, 165114 
(2006). 

25 S. S. Kancharla, B. Kyung, D. Senechal, M. Civelli, M. 
Capone, G. Kotliar, and A.-M. S. Tremblay, Phys. Rev. B 
77, 184516 (2008). 

26 M. Caffarel and W. Krauth, Phys. Rev. Lett. 72, 1545 
(1994). 

27 M. Potthoff, Eur. Phys. J. B 32, 429 (2003). 

28 W. A. Harrison, Electronic Structure and the Properties of 
Solids, The Physics of the Chemical Bond (Freeman, San 
Francisco, 1980). 

29 According to the linear combination of atomic orbitals 
method (Ref. l28h . the hybridization strength between a 
d orbital and a p orbital separated by I is proportional to 
r 7/2 . Thus, taking t ~ 0.5 eV and the NN Cu bond length 
21 ~ 4 A, a rough estimation for a is a ~ 7f/4 ~ 0.9 eV/A 
for the ideal orthorhombic distortion u. 

30 F e i e vs St analysis becomes unreliable at TV < 0.73 because 
the distance function which emphasizes the low-frequency 
regime proportional to 1/|oj„| is not optimal for the over- 
doped Fermi-liquid regime.— The current choice of the 
distance function results in the increase of F e i e with 8t- 
However, it is expected that the phase boundary for fi- 
nite U will asymptote to that for U — with decreasing TV 
when an appropriate distance function for the Fermi-liquid 
regime (Ref. [3lh or a more reliable impurity solver such as 
the continuous-time quantum Monte Carlo (Ref. I32T ) is em- 
ployed. 

31 D. Senechal, Phys. Rev. B 81, 235125 (2010). 

32 E. Gull, A. J. Millis, A. I. Lichtenstein, A. N. Rubtsov, M. 
Troyer, and P. Werner, Rev. Mod. Phys. 83, 349 (2011). 

33 D. Senechal, D. Perez, and M. Pioro-Ladriere, Phys. Rev. 
Lett. 84, 522 (2000). 

34 I. Khavkine, C.-H. Chung, V. Oganesyan, H.-Y. Kee, Phys. 
Rev. B 70, 155110 (2004). 

35 F. Laliberte, J. Chang, N. Doiron-Leyraud, E. Hassinger, 
R. Daou, M. Rondeau, B.J. Ramshaw, R. Liang, D.A. 
Bonn, W.N. Hardy, S. Pyon, T. Takayama, H. Takagi, I. 
Sheikin, L. Malone, C. Proust, K. Behnia, and L. Taillefer, 
Nat. Commun. 2, 432 (2011). 

36 R. S. Markiewicz, J. Lorenzana, G. Seibold, A. Bansil, 



6 



Phys. Rev. B 81, 014509 (2010). 

R. S. Markiewicz, J. Lorenzana, G. Seibold, A. Bansil, 
larXiv: 1207.57151 

H. Yamase and W. Metzner, Phys. Rev. Lett. 108, 186405 
(2012). 

N. Poccia, G. Campi, M. Fratini, A. Ricci, N. L. Saini, and 
A. Bianconi, Phys. Rev. B 84, 100504(R) (2011). 

I. Zeljkovic, E. J. Main, T. L.Williams, M. C. Boyer, K. 
Chatterjee, W. D. Wise, Y. Yin, M. Zech, A. Pivonka, T. 
Kondo, T. Takeuchi, H. Ikuta, J. Wen, Z. Xu, G. D. Gu, 
EW. Hudson, and J. E. Hoffman, Nat. Mater. 11, 585 
(2012). 

A. Migliori, William M. Visscher, S. E. Brown, Z. Fisk, S.- 
W. Cheong, B. Alten, E. T. Ahrens, K. A. Kubat-Martin, 
J. D. Maynard, Y. Huang, D. R. Kirk, K. A. Gillis, H. K. 
Kim, and M. H. W. Chan, Phys. Rev. B 41, 2098 (1990). 
J. Wu, Y. Wang, H. Shen, J. Zhu, Y. Yan, and Z. Zhao, 
Phys. Lett. A 148, 127 (1990). 

D. Fournier, G. Levy, Y. Pennec, J. L. McChesney, A. 
Bostwick, E. Rotenberg, R. Liang, W. N. Hardy, D. A. 
Bonn, I. S. Elfimov, and A. Damascelli, Nat. Phys. 6, 905 
(2010). 

N. Lin, E. Gull, and A. J. Millis, Phys. Rev. B 82, 045104 
(2010). 



45 H. Yamase and R. Zeyher, Phys. Rev. B 83, 115116 (2011). 

46 A. W. Rost, R. S. Perry, J.-F. Mercure, A. P. Mackenzie, 
and S. A. Grigera, Science 325, 1360 (2009). 

47 E. C. Blomberg, M. A. Tanatar, R. M. Fernandes, 
B. Shen, H.-H. Wen, J. Schmalian, and R. Prozorov, 
larXiv:1202.4430l 

48 C. Honerkamp, Phys. Rev. B 72, 115103 (2005). 

49 B. Binz and M. Sigrist, Europhys. Lett. 65, 816 (2004). 

50 S. Raghu, A. Paramekanti, E.-A. Kim, R. A. Borzi, S. A. 
Grigera, A. P. Mackenzie, and S. A. Kivelson, Phys. Rev. 
B 79, 214402 (2009); W.-C. Lee and C. Wu, Phys. Rev. B 
80, 104438 (2009). 

51 C. Fang, H. Yao, W.-F. Tsai, J. P. Hu, and S. A. Kivelson, 
Phys. Rev. B 77, 224509 (2008). 

52 R. M. Fernandes and J. Schmalian, Supercond. Sci. Tech- 
nol. 25, 084005 (2012). 

5,i F. Kriiger, S. Kumar, J. Zaanen, and J. van den Brink, 
Phys. Rev. B 79, 054504 (2009). 

54 A. Cano, M. Civelli, I. Eremin, and I. Paul, Phys. Rev. B 
82, 020408(R) (2010). 

55 S. Liang, G. Alvarez, C. §en, A. Moreo, E. Dagotto, Phys. 
Rev. Lett. 109, 047001 (2012). 



